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Abstract 

We study two simple real analytic uniformly hyperbolic dynamical systems: ex- 
panding maps on the circle 5 1 and hyperbolic maps on the torus T 2 . We show that the 
Ruelle-Pollicott resonances which describe time correlation functions of the chaotic 
dynamics can be obtained as the eigenvalues of a trace class operator in Hilbert space 
L 2 (S 1 ) or I? (T 2 ) respectively. The trace class operator is obtained by conjugation 
of the Ruelle transfer operator in a similar way quantum resonances are obtained in 
open quantum systems. We comment this analogy. 

PACS numbers: 05.45.-a , 05.45.Ac, 05.45.Pq 
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1 Introduction 

In this paper, we study the Ruelle-Pollicott resonances of two simple models which are 
uniformly hyperbolic: expanding maps on the circle S l and hyperbolic maps on the torus 
T 2 . These models are the most simple examples of chaotic dynamical systems, which 
exhibit strong chaotic properties, such as ergodicity, mixing, decay of correlations, central 
limit theorem for observables, etc. , see [14] [9] [3]. Expansivity or hyperbolicity makes every 
trajectory unstable and are therefore important hypothesis responsible for these chaotic 
properties. One of these properties, the "decay of time correlation" is fundamental to 
establish other chaotic properties. Decay of time correlations can be studied by means of 
the spectral analysis of "Ruelle transfer operators", which transport functions on S* 1 or T 2 
according to the dynamics [7, 20]. The simplest Ruelle transfer operator is the "pull back 

operator" or "Koopman operator" defined by: ( Mip ) (x) = f ip (M (x)), where M : S 1 — > S 1 
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is the map (resp. M : T 2 — > T 2 ), and <p G C° (S 1 ). The time correlation of two functions 
(/>,ip E C°(S l ) is defined by C^(t) = (^|M*b) = J sl (j) (x) (<p o M l (x)) dx. The main 
effect of the hyperbolic dynamics is to transform an initial function tp into a function M t ip 
with finer and finer structures, as the time t evolves. In other words, the information on 
the initial function <p is sent towards "microscopic scales", or equivalently at infinity in the 
Fourier space. On the "macroscopic scale" (i.e., if M f ip is tested on a regular function 0), 
there remains only a constant function, i.e., the function 1 times a weight Ushb (v 9 ) where 
I^srb is called the Sinai-Ruelle-Bowen (S.R.B.) measure. The number Hsrb (<f) G C is the 
only information on <p which remains on the macroscopic scale after a long time evolution. 
With this point of view, the decay of time correlation functions (t) is due to the escape 
of the function M l ip towards infinity in the Fourier space, implying a decay of the small 
Fourier components. This suggests to study the decay of correlations using a "window of 
observation" in the Fourier space, centered on small Fourier components (in the unstable 
direction). This is the role of the operator A below. This situation is very similar to open 
quantum systems, where the decay of the quantum wave function in a compact domain 
of space is due to the escape of the wave function towards infinity. In such situations 
people study the decay by a "complex scaling method" which consists in conjugating the 
dynamical operator in such a way that the wave function is toned down at infinity ([6], 
chap. 8). Then, the "quantum resonances" which appear are suitable to describe the decay. 
Here we will follow this general approach. 

We will suppose in both models that the map M is real analytic. We will show that 
the time correlation functions can be obtained from an "effective dynamical operator" R 
obtained from M by a conjugation R = AM A' 1 , where A tones down the high Fourier 
modes in the unstable direction of the dynamics. The main result is to show that R is a 
trace class operator in the Hilbert space L 2 (S l ) (resp. L 2 (T 2 )). The effective long time 
dynamics is obtained by R = AM t A~ 1 with t 6 ff, and the spectral properties of R are 
important for that. The eigenvalues of R are called the "Ruelle-Pollicott resonances". This 
approach, with a conjugation, is the one which is usually followed in the complex scaling 
methods [6]. An alternative but equivalent approach would be to keep the operator M but 

to work with another norm instead of the L 2 one, namely with ((p,ip) A == (Acj) , Aipj ^. 

This last formulation is preferred in [5] [11] [15] [4] . In this approach, the operator A is seen 
as an isomorphism between two Hilbert spaces, and one gets yf), Mij^j = Rip'^j ^ 

with ft = A(j>, i)' = Aip. 

The correlation spectrum for analytic maps has already been studied through trace class 
operators but with different approaches: the case of expanding maps has been studied by 
D. Ruelle in [20]. The case of analytic hyperbolic maps has been studied by H.H. Rugh 
in [21]. Our approach which consists in working in the Fourier space (or more technically 
conjugating M by a pseudo-differential operator) has been already investigated recently for 
hyperbolic diffeomorphisms by V. Baladi and M. Tsujii [4]. In their paper, they consider a 
much broader class of dynamical systems (they do not suppose analyticity) and they show 
quasi-compacity for the transfer operator in a Banach space of distributions. Although our 
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results are more restrictive, we believe that they have their own interest because of their 
technical simplicity (we obtain a trace class operator in a Hilbert space and the proof is 
quite simple). A technical difference between the two approaches appears for example with 
the choice of the operator A. We have to choose an operator which has an exponential 
expression in the Fourier basis, whereas in [4] the operator A has an algebraic dependence 
on the Fourier modes (A is a power of the Laplacian). This difference is crucial to obtain 
some of the results presented in this paper. 

Our techniques do not allows us to treat any hyperbolic maps on the torus, but maps 
which are closed enough to the linear hyperbolic map 1 . More precisely, as the proof will 
show in Section 4.2, the technical assumption is that the unstable and stable tangent 
directions are uniformly contained in constant cones defined by the linear map. This 
restriction of our method, due to the very simple expression of the operator A, saves us on 
the other hand from making a partition of the unity as in [4]. 

With a different approach, results close to those obtained by V. Baladi et al. were ob- 
tained by C. Liverani et al. in [5][11][15]. The connection between dynamical determinants 
and Ruelle resonances is established there in great generality. See [4] for historical remarks 
and comparison between the different approaches. 

The paper is organized as follows. In Section 2, we define the operator R which governs 
the decay of the correlation function. We state the theorems which say that R is a trace 
class operator in both case S 1 and T 2 . These results are proved respectively in Section 
3 and 4. Our method is very convenient for numerical analysis since the truncation of 
the high Fourier modes produces a small error. In Section 3.3, we present numerical 
illustrations of some aspects of the S.R.B. measure and the R.P. resonances. In Section 
2.1.4, we show the equivalence of our approach with a more common approach known as 
"randomly perturbed dynamics" or "noisy models" [8] [5]. We use this equivalence to show 
the (well known) Trace formula in terms of periodic orbits [3]. We mention that in [11], 
a powerful method is developed in the C r case and allows the authors to show spectral 
stability for a wider class of deterministic and random perturbations. 

Acknowledgement. We would like to thank Viviane Baladi, Patrick Bernard, Alain Joye, 
Malik Mezzadri and Stephane Nonnenmacher for discussions related to this work. FF 
gratefully acknowledges partial support by "Agence Nationale de la Recherche" under the 
grant JC05_52556. 

2 Statement of the results 

2.1 Expanding map on the circle 

Let M : E -* E be the map defined by 2 

M (x) = 2x + / (z) 
1 We discuss some possible extension in the conclusion 

2 More generally we could have supposed that M (x) = dx + f (x) with d £ N, d > 2. This does not 
change the results we obtain. 
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where / is real analytic and periodic: f(x + l) = f (x) ,Wx G R. We suppose moreover 
that 

def . / df\ 

/min = ^ XT > "! 



<7.r 



so that M' x = 2 + -£ > 1 (M is called strictly expanding). A simple example used later for 
numerical illustrations is 

/(x) = ^-sin(27rx), |5|<1 (1) 

For all n e N and all x G R, M (x + n) = M (x) modi, hence M defines an expanding 
map on the circle S 1 = R/Z (also denoted by M). M is not invertible, but it is the 
simplest model exhibiting chaotic dynamics, see [14]: expansivity is responsible for the 
high sensitivity to initial conditions, mixing and positive entropy of the dynamics. 

The pull-back operator M on L 2 (S l ) is the (non unitary) bounded operator defined 

by 3 

(Mcp\ (x) = <p (Mi) (2) 

For any n G Z, we denote 4 by \n) G L 2 (S 1 ) the Fourier mode (p n (x) = exp (i27rnx), 
and define the operator A by 

A\n) = exp (—a |n|) |n), n G Z, with a > 0. (3) 

A is diagonal in the Fourier basis. The image Ca = A(L 2 {S 1 )) is a set of very regular 
functions (analytic in a complex neighborhood of S 1 of radius a). The operator A is used 
to define 5 the operator: 

R = AM A^ 1 (4) 

with domain C^, dense in L 2 (S 1 ). 

In this paper, we will prove the following theorem: 



Theorem 1. There exists a > entering in eq.(S), such that matrix elements of R 



decrease exponentially: 



(n'\R\n) < exp(— c(\n'\ + |n|)), with c > 0. In particular R 



extends to a trace class operator in Hilbert space L 2 (S 1 ] 



3 We can also consider a more general class of operators, called "Ruelle transfer operator": (M g ip\ (x) d = 

e 9 ^ <p (M x) , with complex valued function g. The results obtained in this paper extend to this case 
provided g is real analytic. 

4 Throughout the paper, we use Dirac notations for vectors in Hilbert space H — L? (S 1 ): tp G H is 
written \<p). Its dual metric is written (<^|. A scalar product is written (<j>\tp) = J sl <j> (x)<p (x) dx. If M 

dcf A 

is an operator, we write (4>\M\ip) = (<j>\Mip) = (M*cf)\(p), (where M* is the adjoint operator). Finally, 

\m<p\ = \4)®(<Pi 

5 Equivalently A can be seen as a change of norm on 1? (S 1 ). 
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As the proof will show, the result holds for any a £ ]0, ao], with some a Q > 0, and c 
depends on a. 

For general results on trace class operators, see [18, 19], or [10] chap. 4. The eigenvalues 
of R are called the Ruelle-Pollicott resonances (R.P.) of the pull back operator M. 

2.1.1 Dynamical correlation functions 

The operator R is useful to study time-correlation functions. Indeed, if \<p) £ Ca is a 
regular "test function", and \(p) £ L 2 (S 1 ), then for t £ N, C^ )tp (t) = {^)\M t \ip) can be 
expressed using R as 

(t) d = f (0|M'|^> = ((^li- 1 ) # (A\ V )) (5) 

In "physical terms" it means that the operator R is a nice effective operator to express the 
dynamics of M in L 2 (S 1 ) provided it is tested on the regular function space Ca- 

2.1.2 Finite rank approximation 

There is a direct consequence of Theorem 1 which is useful for the numerical computation 
of Ruelle-Pollicott resonances. Let M N be the matrix of the operator M expressed in the 
Fourier basis and truncated to the N first Fourier modes (i.e., |n| < N , the matrix M N 
has thus size (2N + 1) x (2N + 1)). 

Corollary 2. The eigenvalues of M N converge towards the Ruelle-Pollicott resonances, 
when N —> oo. 

Proof. If Rn is the matrix of the operator R restricted to Hn = Span {\n), \n\ < N} = 
C 2N+1 , then the eigenvalues of Rn converge to the R.P. resonances as iV — > oo, because 
Rn converges to R in operator norm. But in Hn = C Ar+1 , Rn is conjugate to M N by the 
invertible and diagonal matrix A N = Diag (exp (—a \n\) , n = —N — > A r ), R N and M N have 
the same spectrum. □ 



2.1.3 Exponential concentration of R.P. resonances near zero 

Since R is a compact operator (and moreover a trace class operator), its eigenvalues con- 
verge to zero. From Theorem 1, matrix elements (n'\R\n) decrease exponentially fast for 
large \n\ , \n'\. A quite direct consequence of this is the exponential concentration of the 
eigenvalues of R near zero: 

Theorem 3. Let A.; £ C, i = 0,1,... be the non zero eigenvalues of R (the Ruelle- 
Pollicott resonances), such that \\%+\\ < \\\, and counting multiplicity. Let li = log|Aj| 

l( l+e -c ") 

and C\ = , _J . For any i > 0, 

(1 e ) 

Z*< --z + logd (6) 
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The constant c is given in Theorem 1 and the proof of this theorem is postponed to 
Section 3.3. (The best estimate is for the largest possible value of c, which is neither very 
explicit in the proof, nor very sharp). 

2.1.4 Relation with randomly perturbed operators or noisy models 

A different approach to obtain Ruelle-Pollicott resonances of transfer operators is to add 
a small "noise", or "random perturbation" to the dynamical operator M at each step of 
evolution. These models are often used for theoretical or numerical calculations. In [8] 
Baladi and Young show that for expanding maps, the randomly perturbed operator is 
compact, and that its eigenvalues are the Ruelle-Pollicott resonances in the limit when 
the perturbation vanishes. A similar result is shown by Blank Keller and Liverani in [5] 
for Anosov maps. See also [17]. In this section we consider such "noisy operators" and 
show with a very simple proof that they have the same eigenvalues as R (i.e., the R.P. 
resonances) when the perturbation vanishes. The same result holds for Anosov maps on 
T 2 considered in the next section. 

Let A = — be the Laplacian operator on S 1 , and for e > 0, let 

£> £ dcf e -eA/(27r) 2 

be the heat operator. This operator is diagonal in the Fourier basis: 

V £ \n) = e- £n2 \n) (7) 

The main effect of T> £ is to truncate the high Fourier components. In the "real space" 
i 6 S 1 , the operator T> £ convolutes with a Gaussian distribution of size ~ y/e, so T> £ as 
the same effect Gaussian noise. 

Let us define the "noisy perturbed operator" by 

M £ = MV £ (8) 

M £ is a Trace class operator because it is the product of T> £ which is Trace class with M 
which is bounded (cf [18] p.207). 

Theorem 4. The eigenvalues of the noisy perturbed operator M £ converge to the Ruelle- 
Pollicott resonances, when £ — > 0. 

Proof. Let us define the operator 

R £ = RV £ = AMA~ 1 V £ = AMV £ A~ 1 = AM £ A~ x 

where we have used the fact that A and T> £ commute. Let R £t N (resp. M £> n) the matrix 
of the operator R £ (resp. M £ ) expressed in the Fourier basis and truncated to the first A^ 
Fourier modes. We have R £i n = A^^A' 1 so the matrices R £ ^ and M £tN have the same 
spectrum. But R £< n converges to R £ in operator norm for A^ — >• oo (resp. M £>N — > M £ for 
A^ — » oo). We deduce that R £ and M £ have the same spectrum. Now R £ converges to R 
in operator norm for e — >• 0, which gives that eigenvalues of R £ converge to eigenvalues of 
R. □ 
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2.1.5 Trace formula 

An important feature of Ruelle transfert operators is the existence of exact trace formulas 
in terms of periodic orbits ([3] p. 100). We just recall here this trace formula for the 
operator R defined by Eq.(4). 

Proposition 5. For any t > 1, 

1 



ILLM* - II 



x£Fm(M*) 

where Fix(M t ) denotes the set of fixed points of M l , and D x M t (x) = (x). 

Proof. We consider first the trace of the operator M t>£ = f M l t) £ , (similar to Eq.(8)). From 
([10] Th. 8.1 p. 70), 

Tr = dx(x\M t Jx) = J dx5 £ (M* (x) - x) 

where (x'\M tt£ \x) denotes the Schwartz kernel of the operator M t ,e on L? (S 1 ) and where 
S £ = D £ S denotes the regularized Dirac distribution at x = (i.e. 5 £ is a periodic Gaussian 
function with width ~ s/e). With the ((2* — l)-valued) change of variable y = M l (x) —x, 
we obtain 

so Tr (Mt^ — > Y^xeFixiAP) \d x m^i\ for e ~* °- 0n the otner hand, with R t:£ = f £*£> £ , we 

show (as in the proof of Theorem 4), that Tr (M t , e ) = Tr \RtJ\ > and that Tr ~ ~> 
0, for e — > 0. This gives the result. □ 

2.2 Hyperbolic map on the torus 

With the same approach, one can study a non linear hyperbolic map on the torus (expressed 
as a linear map with a small perturbation). Let M G SL (2, Z) be an hyperbolic matrix 
(i.e., TrM > 2) and / : M 2 — > IR 2 be a real analytic periodic function: 

/ (x + n) = f (x) , Va; G IR 2 , G Z 2 

Let M : IR 2 -> IR 2 be defined to be M perturbed by /: 

M (x) = M (x) + / (x) 

Then M (x + n) = M (x) [Z 2 ] hence M induces a map on T 2 also denoted by M. Structural 
stability asserts that the map M on T 2 is Anosov (uniformly hyperbolic) whenever 
is small enough ([2] p.l22)([14] p. 89) 6 . 

6 In our case / is supposed to be real analytic. So ||/||pi is controlled by ||/||co. 
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The pull-back operator M on L 2 (T 2 ) is the bounded operator defined by 

[M^j (x) = if (Mx) . (9) 

Note that the operator M is not unitary except if M preserves the area on T 2 . 

/ 2 1 \ 

Example: Choose M = I - - J and 

/ (x) = (o, ^ sin (2vr (2 Xl + x 2 )fj . (10) 

In this example, M preserves 7 area dx\dx 2 . 

For each n = (rii,ri2) G Z 2 , denote by \n) G L 2 (T 2 ) the Fourier mode (p n (x) = 
exp (i2n (n.x)). Let u, s G M 2 be unstable/stable eigenvectors of the transposed matrix 
M l , i.e., M 1 m = e x °u and M l s = e" A °s, with A > 0. A vector v = (v x , v 2 ) G M 2 is written 
v = (v u , v s ) with respect to the basis (u, s). In particular n G Z 2 is mapped to n = (n u , n s ). 
Define the operator A by 

A\n) == exp (— a \n u \ + a |n s |) |n), nGZ 2 , with a > 0. (11) 

It is diagonal in the Fourier basis. A is defined on a domain Da == Dom = 

{\<P) = J2nfn\n), s.t. Y,nWn\ 2 e 2a \ n ^- 2a \ n -\ < oo, EJ^ 2 < oo} dense in L 2 (T 2 ), and 
consists of functions with exponentially fast decreasing Fourier components (i.e very regu- 
lar) in the stable direction. Similarly, 



„ def 
U A = 



Dom (i- 1 ) = j|0> = s.t. J2\<Pnf e~ 2a ^ +2a ^ < oo, ^ |0„| 2 < oo 

^ n n n 

C L 2 (T 2 ) 

consists of functions with exponentially fast decreasing Fourier components in the unstable 
direction. One checks that Ca = A (Da), Da = Ar 1 (Ca). 
Define 

R = AMAT X (12) 
on a suitable domain included in Ca (one has Dom (^Rj = Ca if M (Da) C Da). 



Theorem 6. For a C l -small enough perturbation f, i.e., \\f\\ c i < s with e > 0, there 
exists a > such that R is defined on the domain Dom (ji^j = Ca and its matrix ele- 



ments decrease exponentially: 



(n'\R\n) < exp ( — c (1^1 + |n' 2 | + |ni| + I^D), with c > 0. 



Therefore, R extends to a trace class operator in Hilbert space L 2 (T 2 



7 Because in this example, M can be written as M = M\Mq, where Mi is an Hamiltonian time 1 flow, 
generated by Hamiltonian function Hi (xi,X2) = t^tz cos (2-kxi) on T 2 . 
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As the proof will show, the result holds for any a £ ]0, ao], with some a Q > 0, and c 
depends on a. The operator R is useful to study time-correlation functions: if \<f>) £ Ca is 
a regular "test function", and \tp) £ Da, then 

it) ¥ <0|M» = # (i|^>) 

The different corollaries and applications we mentionned for expanding maps, work as 
well for hyperbolic maps on the torus (with only minor modifications): (1) finite rank 
approximation, (2) exponential concentration of R.P. resonances near zero, (3) relation 
with randomly perturbed operators, (4) Trace formula in terms of periodic orbits. 



2.2.1 Exponential concentration of R.P. resonances near zero 

Theorem 7. Let X n £ C, n = 0, 1, . . . be the non zero eigenvalues of R (the Ruelle-P ollicott 
resonances), such that |A„ +1 | < |A n |, and counting multiplicity. Let l n = log|A n |. There 
exists C\ > such that for any n > 0, 

the constant c is given in Theorem 1, and the proof of the theorem is postponed to 
Section 4.2.1. The main difference with (6) is the appearance of n 1 ^ 2 , where the power is 
1/d with d = dim (T 2 ) = 2. Let us remark that such an asymptotic behaviour of eigenvalues 
is also met in quantum mechanics in the spectrum (E n ) n of the Harmonic Oscillator on 
R. d . From the semi-classical Weyl law E n ~ Csten 1 ^ (see [13] chap. 16). This suggests 
that (13) could be obtained or interpreted within a semi-classical approach, with a Weyl 
asymptotic. 



2.2.2 Relation with randomly perturbed operators 

The analysis made in section 2.1.4, can be repeated with no change except for the definition 
of the heat operator: 

f) £ \n) = e -< n * +n *)\n) 
which is used to define the randomly perturbed operator: 

M £ = MV £ 

We have 

Theorem 8. The eigenvalues of the noisy perturbed operator M e converge towards the 
Ruelle-P ollicott resonances, as £ —>■ 0. 

(The same proof as in section 2.1.4.) 
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2.2.3 Trace formula 

We have the following Trace formula for R l in terms of periodic orbits. 
Proposition 9. For any t > 1, 

1 

I det ( D ,. 



\det{D x M l - Id)\ 



where Fix(M t ) denotes the set of fixed points of M t , and D x M l (x) is the differential at 
point iGT 2 . 

The proof follows the same lines as for Proposition 5. 

3 Expanding map on the circle 

In this section, we prove Theorem 1. 

3.1 Matrix elements of the operator M 

Denote \n) the Fourier mode ip n (x) = exp (i2irnx), with n G Z, x G S 1 . The set (\n)) neZ 
form an orthonormal basis of L 2 (S 1 ) and matrix elements of M in this basis are explicitly 
given by 

(n'\M\n) = / dx exp (i2n ((2n - n) x + nf (x))) (14) 



3.1.1 Remarks 

• For a vanishing perturbation / = 0, then 

(n'\M \n) = 5 2n=n , (15) 

i.e., in the plane (n',n), matrix elements are zero except on the "line" n' = 2n. For 
a non zero perturbation /, we will show that matrix elements are very small outside 
a cone containing this line. 

• Since / is real, we have the symmetry 



(-n'\M\ - n) = (n'\M\n) 

and if n — 0, we have 

(n'\M\0) = 5 n , =0 
Therefore we have only to study matrix elements with n > 0. 
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3.1.2 Localization property of matrix elements 

For simplicity of the presentation, we borrow notations from semi-classical analysis (see 
[16]). For n > 0, let us make the change of variables (n, n') <^> (h,p), with 

h = — , p = — On! — 2n) = nh — 2 (16) 
n n 

and define K=h/ 02n). A matrix element can be written as the oscillating integral: 



In Op) = f (n'\M\n) = I dx exp (i (/ (x) — px) /K) 
Jo 

From the "non stationary phase theorem" below, this integral is "very small" for values 
of p outside the interval [ / m in , / max ], with 



H dcf . df , dcf df 



Theorem 10. "Non stationary phase". Assume that f (x) is a periodic function and 
can be continued analytically in some strip \Im{x)\ < Y. Assume pj 02itK) G Z, and 
In (p) = J Q dxexp (i (/ (x) — px) /K). For any e > 0, there exists a constant C > 0, such 
that for any p, K > 0, 

\h (p)\ < min (l, e" c ( / — - p " £ ) /a , e - cr (P-/m-- e )/ R ) (17) 





£ 










1^ 
















1 




- 

P 



X/ U £1 

J rain J max 

Figure 1: Upper bounds for the function \In{p)\ 



Proof. Write z = x+iy and f (z) — a (z)+ib (z), with a, 6 real valued functions. Analyticity 
of / gives d y b = d x a. For y = and x G [0,1], one has b (x, 0) = and (c^fe) (x, 0) = 
(3e°)j,=o > /min- Therefore Ve > 0, 3y > y < Y such that b (x, y ) > (/ min - e) y for 
all x. So for 2 = x + iyo, 

|exp (i (/ (z) - pz) /fr)| = exp (- (b (x, y Q ) - py ) /h) < exp (- (f min - e - p) y /h) . 

Since / is analytic, we can deform the integration path to z = x + iyo, x — —> 1, yo fixed. 
This gives 

\h (p) I < exp (- (/ min -e-p) y / h) . 
A similar argument for the integral In Ojp) = Jq dxexp (i (— / (x) + px) /K) provides 

14 (p) I < exp (- (p - f max - e) y'Jh) 
We finally choose C = min (y , y' Q ). □ 
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3.1.3 Remarks 

• For p e [/mm) /max]) the stationary phase formula gives the asymptotic value of In (p) 
when h — > 0. It says that the asymptotic value of the matrix element (n'\M\n) = 
In (p) of operator M depends only on the point x' in the integral, such that P = (%'), 
equivalently n! = (2 + £ (x')) n. We comment now on a semi-classical interpretation 
of this result. The map M _1 acting in S 1 is two valued. Let M~ l denotes its 
lifted action on the cotangent space T^S 1 . If x = M(x') = 2x' + f (x'), then = 
(2 + £ 0')) £. So, if (x, k) denote coordinates on T*^ 1 , and (x', k') = M~ l (x, k), 
then x = M (x') and fc' = (2 + £ (x')) fc. So in a sense which need to be precised, 
M acting in L 2 (5* 1 ) is a "semi-classical quantization" of the map M _1 acting in the 
symplectic space T*^ 1 . This "semi-classical aspect" of hyperbolic dynamics will be 
investigated in a futur work. 

• If we come back to variables n,n', this last theorem shows that matrix elements 
(n'\M\n) in the plane (n',n) are exponentially small outside the cone defined by (for 
n > 0): 

(f^ Q + 2-e)n<n' <(f^ ax + 2 + e)n 

Expansivity hypothesis of M gives f min + 2 > 1, so for e small enough this cone does 
not contain the diagonal n' = n. Cf Figure 2. 




n 



Figure 2: Matrix elements (n'\M\n) are small outside the grey cone, for |n| — ► oo. For 
/ = 0, this cone reduces to the line n' = 2n. These matrix elements can be interpreted 
as transition amplitudes for the dynamics n — > n' . For long times (many iterations), the 
dynamics goes to infinity in Fourier space on the sector \n'\ > \n\, and tends to small 
Fourier components on the sector \n'\ < \n\ (see the black arrows). Important matrix 
elements of M are localized on the sector \n'\ > \n\ only, and thus generate an escape 
towards infinity in Fourier space, responsible for chaos, as discussed in the introduction. 
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With the example (1), we can explicitly express matrix elements in term of Bessel 
functions of the first kind (see [1] 9.1.21 p. 360): 



(n / \M\n} = (-lf n - n ' ) J {2n „ n , ) (8n) 



This gives 



\h(p)\ 



J(- P /h) 



Asymptotic results for Bessel functions (see [1] 9.3.1 and 9.3.2 p. 365) give 

log|ift(p)| ~ — ^ ^log ^ ' for /i fixed, and p ^ oo 

p p 
~ — — (a — tanh a) , with cosh a = - > 1 fixed, and h — > 
h o 

which are sharper than the upper bound in (17). 



3.2 Proof of Theorem 1 

3.2.1 Idea of the proof and remarks 

In the linear case, with a vanishing perturbation / = 0, the result is obvious. Indeed 
from eq.(15), matrix elements of R lie on the line n' = 2n and decrease like (n'\R\n) = 



<L/ = 



a(|n|-|n'|) 



S n '= 



:2ne 



-a|n| 



(the spectrum is then a [Rj = {1} U {0}, with 1 as a 



simple eigenvalue). Remark that choosing the operator A with the algebraic form A\n) 
l^pr|n), n G Z would not give decreasing matrix elements: (n'\R\n) 



* , |n|a 

u n—2n | n /| a 



^n'=2n^r- On the other hand, the choice A\n) = e ' n ' |n) with < a < 1 or A\n) = 
e -aio g 2 (i+|n|)| n ^ wou i(j be suitable as well. 

In the non linear case, with / ^ 0, we have shown that (n'\M\n) decreases fast 

outside a cone in the (n',n) plane. The conjugation with A gives the multiplicative factor 
e+ a d n H n l) s which decreases in the sector \n'\ > |n|, but increases in the sector \n'\ < \n\. 
The expansivity hypothesis insures that the cone is strictly included in the first sector. 



Moreover, provided a > is small enough, the decrease of 



(n'\M\n) 



dominates the increase 



of e + a (M l n 'l) j n the sector Ira'l < Iral. At final result we obtain that 



exponentially fast for \n\ 



{n'\R\n) 



decreases 



\n 



oo. 



3.2.2 Exponential decrease of matrix elements 

From (3) and (4) it follows that for any n > 0: 

(n'\R\n) = e +a{n - ln ' l) (n'\M\n). 
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Instead of n',n G Z, we prefer to use "re-normalized" indices for n > 0, defined by 

v' = nh = p + 2 



?7 



where were already defined in (16). 

We can write e + a ( n -\ n '\) = e^ A(v '\ with Aiv 1 ) = f a(l — |z/|), and Equation (17) gives 
the upper bound: 



(ra'|M|n) 



< c 



B (z/) d ^ f min (0, 2ttC (z/ - 6 min ) , 2vrC (6 max - z/)) 



with 



2 + /min £ > 



^ + /max + £ - 



Because M is expanding, we can choose e > such that min x (M^) = 2 + / min > 1 + e, 
hence b m > 1. For large |z/|, the functions A (z/) and B (i/) and have respective slope 
a and 27rC. Choosing a < 27rC implies that the maximum of F (z/) = A (z/) + B (z/) is 
reached for v' — b m . See Figure 3. 



-c(i + M) 




F = A + B 



Figure 3: Representation of functions F (z/) = A (z/) + 5 (z/) for the upper bounds 
f\R\n) <exp(|F(z/)) < exp (-f (1 + |z/|)) . 



n 



We need an upper bound 



(n'\R\n) < exp (— c(\n\ + |n'|)) = exp (— | (1 + |z/|)), there- 
fore we look now for a constant c > such that F (z/) = A(z/) + B < — c (1 + |z/|), 



< a. 



for any z/ G K. This requires c < a and F (& mm ) < — c(l + b m - m ) ^ c < a 
Consequently we choose c < a S in ~! , and this proves the exponential estimates of Theorem 

min ~> 



6min + l 



1. 



3.2.3 Trace class operator 

First from (n'\R\n) < exp (— c{\n\ + \n'\)), R is a Hilbert-Schmidt operator, therefore 

bounded. Its domain Ca is dense in L 2 (S 1 ). From a classical result, -R extend in a unique 
way to a bounded operator in L 2 (S 1 ). Now let B be the operator diagonal in the Fourier 
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basis, defined by B\n) = -nbt\n), with a > 1/2. B is a Hilbert-Schmidt operator and 

C = f RB ^ is also an Hilbert-Schmidt operator, so R = CB being a product of two 
Hilbert Schmidt operators is a trace class operator (cf. [10], Lemma 7.2, p. 67). 



3.3 Exponential accumulation of R.P. resonances near zero 

In this section, we prove Theorem 3. 



[n'\R\n) 



< 



Let R be the trace class operator obtained in Theorem 1, with the estimation 

exp (— c(\n\ + \n'\)) on its matrix elements, with c > 0. We first deduce an estimation on 
the singular values of R: 

Lemma 1. Let j = 0,1, ... be the non zero singular values of R (namely the eigenvalues 

of the self-adjoint operator R), such that Hj+\ < [ij, repeated as many times as the 
value of their multiplicity. Then 

lij < C ie - Cj/2 (18) 

with Ci = ( i_ e - c) 2 • 

Proof. We borrow an argument from [12] (in the proof of Prop. 3.2). From the min-max 
theorem, 



mm max 

VcL 2 ^ 1 ), codimV=j vGV, \\v\\=1 



Rv 



Consider the Fourier basis \n),n € Z and Vi = span (|n)), n | >z , hence codimVJ = 21 + 1. If 
\v) G Vj, we compute 



R\v) 



\n')(n'\R\n)(n\v) 

n'eZ,\n\>l 



— | ( n '\R\ n ){ n \ v ) 

n'eZ,\n\>l 



-c\n\ 



E 



with Sj 



dcf 



— IMI ex P (~ c (l n l + \ n '\)) = \\ v \\ I e ~ " 

n'eZ,\n\>l \H>1 

= \\v\\ 2S l+1 (1 + 25i) = ||v|| e~< l+1) 2S (1 + 25i) 



■c\n'\ 



l-e~ c 



e~ cj S . We deduce that // 2l +i < e~ c ^ l+1 hS (1 + 2^; 



hence for j odd,^ < e - cj/2 e- c/2 2S (1 + 2Si) < C^' 2 with d = 2,So (1 + 2S X ) = 
2 (1 + e^ c ) / (1 - e" c ) 2 . For j even, ^ < ^ < e~ cj/2 2S (1 + 25i) = C ie - cj/2 . □ 

There is a fundamental relation between eigenvalues (Aj) . of R (sorted such that 
< \Xj\, and repeated as many times as the value of their multiplicity) and singular 
values (fjij). (cf. [10], Th. 3.1 p.52): 



j=0 j=0 



n > 



(19) 
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For non zero eigenvalues, define 

lj = log | Xj | , rrij = log fij 

(These sequences tend to — oo as j — »■ oo). Then the above inequality reads YTj=oh — 
Y^j=o m j- Eq. (18) gives rrij < logCi — cj/2. We deduce that 5Z™ =0 lj < (n + 1) logCi — 

| n(n+l) ^ ] lence 

1 " C 

> L < logCi n 

But l n < lj for j < n, so l n < r n+x \ Y^j=o h — l°S^i — f n i which proves Theorem 3. 

3.4 Numerical illustrations: Sinai- Ruelle-Bowen Measure and Ruelle- 
Pollicott resonances 

In order to illustrate the previous result, we discuss here some well-known aspects of the 
S.R.B. measure and Ruelle-Pollicott resonances of Example (1), obtained by numerical 
diagonalization of operator R (in the Fourier basis). 

3.4.1 The Sinai- Ruelle-Bowen measure 

The zero Fourier mode (constant function) \vq) — \n — 0) is an eigenvector of M (and thus 
R) with eigenvalue Ao = 1. It is known that expanding maps such as eq.(2) are mixing 
[14], which implies as we will see that A = 1 is an isolated eigenvalue of multiplicity 1, 
and all other eigenvalues of R are |Aj| < 1, i = 1,2.... Let |iw ) G L 2 (S 1 ) be the dual 
eigenvector, i.e., (v \wq) = 1 and R*\w ) = \w ) <S=> (w \R = (w \. So in operator norm, 

R f = |^ >HI + C(|Air), | Ax | < 1. 

If \(p) G L 2 (S* 1 ), and |0) G Ca, eq.(5) gives exponential decay of correlation for large t: 

C^(t) = ( ( p\M t \ip) = (0|i-> o )K|i|^) + 0(|Ai|*), |Ax|<l 

where \hsrb) = ^l^o) G Ca is called the S.R.B. measure. Its density is a regular (real 
analytic function) on S 1 , cf Figure 4. The physical meaning of the last equation is that for 
large t the function M l \ip) behaves (as seen from test functions, i.e., from a macroscopic 
point of view) like the constant function \vq) times ([J,srb\<p)- This is mixing property. An 
other interpretation of [isrb is that for almost all xq G S 1 , 

1 T 

(^srb\ = lim -yV<Wx | 

1 — >OQ I ' * 

t=l 

(equality of measures) where 5 X is the Dirac measure at x G S 1 . The right hand side is 
called the "physical measure", because it is constructed from a typical trajectory ([7] p. 640, 
[3] p.73), cf Figure 4. 
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Figure 4: S.R.B. measure: the solid line is the density Psrb (#) computed from a numerical 
diagonalization of R (in the Fourier basis), for the perturbation given in (1), with 5 = 0.4. 
The histogram is constructed from a trajectory of length T = 10 7 starting from a random 
initial point Xq. 

3.4.2 The Ruelle-Pollicott resonances 

Suppose for simplicity that the first N eigenvalues of R are simple, R\Vj) = \i\vi), % = — > 
(N — 1), with Ao = 1, |Aj| < 1, |Aj+i| < |Aj|, and |Ajv| < | Ajv— i | - Let us write \wi) the dual 
vectors ,i.e., (wi\vj) = Sij and R*\wi) = \i\wi) (wi\R = \i(wi\. Then 

JV-l 

ft = \v )(w \ + KfaXwl + o (\\ N f) , \X N \ < 1 

i=l 

shows that the Ruelle-Pollicott resonances Aj govern the asymptotic behaviour of the cor- 
relation functions (5) and the convergence towards equilibrium: 

N-l 

(t) = (0|MV) = ((P\v )(fisRB\v) + E X'MA-^iw^) + O (\\ N \ l ) 

i=l 

Xote that in this last expression \vi), \wi) G L 2 (S* 1 ), hence A\wi) G Ca is a regular function, 
but A _1 |ui) may not belong to L 2 (S 1 ). We have to interpret A" 1 )^) as a linear form on 
the space Ca- Vectors \vi), \wi) depend on the choice of operator A, but eigenvalues Aj, 
and distributions A\wi) , A~ x \vi) do not. 

Let us write \ = pie ld \ pi > 0, hence logAj = log(pj) + Figure (5) shows the 
9 first Ruelle-Pollicott resonances log(Aj), obtained by a numerical diagonalization of R 
in the Fourier basis. In particular, we remark that there are two symmetric clusters of 
eigenvalues. A 4 and A 5 are very close from each other: logp 4 = —5.271, # 4 = 1.101, and 
logps = —5.285, #4 = 1.119. We have no explanation for this. In semi-classical analysis 
such clusters occur because of the "tunnelling effect". It would be nice to find such a semi- 
classical interpretation here. Note that (3) predicts that the values log (pA tend at least 
linearly to — oo, as i —> oo. 
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Figure 5: The first nine Ruelle-Pollicott resonances Aj = pit ldi in log scale, in example 
Eq.(l), with 5 = 0.4. We remark two clusters of nearby eigenvalues. 



4 Hyperbolic map on the torus 

We follow essentially the same lines as in the case of an expanding map on the circle, in 
order to prove Theorem 6. 

4.1 Matrix elements of the operator M 

Consider the operator M defined in Eq.(9). Let \n) denote the Fourier mode on T 2 = 
M 2 /Z 2 , defined by ip n (x) = exp {ilim ■ x), with n G Z 2 , x G T 2 . The set (\n)) neZ2 forms a 
orthonormal basis of L 2 (T 2 ), and matrix elements of M in this basis are explicitly given 
by: 

(n'\M\n) = I dx exp (— i2^n' ■ x) exp (i27ra • (Mo (x) + f (x))) 
Jt 2 




dx exp (i2-K ((Mq (n) - n') ■ x + n ■ f (x))) 



with transposed matrix Mq. 



4.1.1 Remarks 

• For a vanishing perturbation / = 0, then 

(n'\M \n) = 6 n , =M t n (20) 

For a non vanishing perturbation /, we will show now that in the plane n' = (n^, n' 2 ) 
matrix elements are very small outside some domain surrounding the point n' = M^n. 



19 



Since / is real, we have the symmetry 



(-n'\M\ -n) = (n'\M\n) 

and if n = 0, we have 

(n'\M\0) = S nl=0 . 

• With the example given by Eq.(10), one can compute explicitly the matrix elements 
in terms of the Bessel functions of the first kind (cf [1] 9.1.21 p. 360): 

<n'|M|n> = (-l) Nl 5 N2=0 J Nl (Sn 2 ), withiV= (n - (M^) 1 („')) = { N*=^+n[-2n> 2 

4.1.2 Localization property of the matrix elements 

Let us make the following change of variables (n, n') <^> (h, u,p), for 

v = - — - E S , h = — >0, p = h(n'-M* (n)) E M 2 , (21) 
\n\ \n\ 

with \n\ = \Jn\ + n 2 and S 1 the unit circle in Fourier space M 2 . Define % — hj (2ir). Any 
matrix element can be written as the oscillating integral 

In,v (p) = (n'\M\n) = dx exp (i (v ■ f (x) — p ■ x) /K) . 



Theorem 11. "Non stationary phase". Let f : T 2 —> M, 2 be real analytic, H > 0, v E 
U (1), and pj (2Trh) E Z 2 . Then v • (D x f) E IR 2 . Consider the compact domain £ = 
{v ■ (D x f) s.t. x E T 2 ,u E S 1 } C IR 2 which contains 0. We denote by \f[ min , /( ma J and 
[f2,miM f2,max\ ^ e projections of £ the axis p\ and P2 respectively . See Figure 6. For any 
e > 0, there exists C > 0, such that for any p = (pi,p 2 ) with p\ < f[ min — e, any fx > 
and any v E S 1 , one has 

\h,u(p)\ < e- c ( f ^- pi - £ ) /h (22) 

Similarly we have exponential upper bounds for the other three half-planes p\ > f[ max + e, 
Pi < f 2 min — e and p 2 > f 2 max + £ with the same constant C . Moreover we have always 
the general bound \I% jV (p)| < 1. 

Proof. Let e > and write (p) = dx 2 e-'W'ft with Xl (x 2 ) = J* dx 1 e K^f(^2)- P ixi)/n_ 
For v, x 2 fixed, let / (xi) = v ■ f (xi,x 2 ). Then it follows from (17) that, 

ft(x 2 )| < e - Cl ( f ~'™- pi - £ ) /h 

where C\ and /^ in depend on v and x 2 . Let f[ min = min,,^ (^-^g^—^j an d C = 

mm v>X2 (Ci) > 0. Then, for Pl < (f[ Mn - e), one has |/ v (p)\ < ft (s 2 )| < e^'W^" 6 )/*. 
Similarly, we define f[ max , f 2 min and f 2 max and obtain similar estimates. □ 
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P2 




Figure 6: Picture of the domain 8 



Remark 

• We have shown that the integral (p) is "small" outside a rectangle containing p = 
in the plane R 2 . This rectangle shrinks to when the perturbation / is C\ small. 
Coming back to variables (n',n), this means that for n fixed, the matrix elements of 
(n'\M\n) are "small" except in a domain surrounding the point n' = MqU. 

4.2 Proof of Theorem 6 

Instead of the variables n, n' e Z 2 , we prefer to use for n ^ 0, 

\n\ \n\ \n\ 

From (21), we have p = v' — M\v. From (11) and (12), we have for 

R n ^ n = (An'lMlA^n) = exp Q^(^)J (n'\M\n) 

with 

A{u') =a(K| - \u s \ - K| + \u' s \), 

where n, h, v are considered as fixed in the discussion. Note that we do not use yet the 
notation (n'\R\n), but R n i ;n instead, because we don't know yet if \n) belongs to the domain 
of R. This will be proved a few lines below. 

Then Equation (22) gives the upper bound: (n'\M\n) < exp (j^B (z/)) , where the 

function B (u') is equal to on a rectangle domain denoted by Zg(z/), containing the 
point MqU. The function B (z/) decreases linearly outside this rectangle, with a slope 2nC 
which does not depend on v and h, see Figure 7. The size of the domain Z B (u) goes to 0, 
whenever — > 0. 
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We deduce that \R n ^ n \ < exp (±F (z/)), with F (z/) = A (z/) +B(v'). The function 
A (is') is zero for v' = z/, and it is negative and decreases with a constant slope a on a 
domain denoted by Za(v), see Figure 7. At the point v 1 = M^(u) = [e x °u u , e~ x °u s ) , 
the value of A (V) = -a (e A ° - l) + \u 8 \ (l - e~ A °)) < A < is strictly negative, 
uniformly with respect to v G S 1 . Therefore, the domain Zb (v) is strictly included in 
Za (y) if the perturbation / is small enough in C 1 norm. 



/ / Zone Z B {y) 




Level sets of B < 



Level sets of A Zone Z A (u) where A < 



Figure 7: Schematic representation of level sets of functions A (z/) and i? (z/), with respect 
to the frame of unstable/stable directions. Here \v\ = 1, and v is sent to Mq [y] by the 
dynamics. 



If we choose a such that A increases slower than B decreases (i.e., a < c'2-nC where 
d > 0) then there exists c > such that 

F(u') < -c(\u x \ + \u 2 \ + \u[\ + \v' 2 \) 

This gives |-R n ' jn | < exp (— c(|n-i| + \n 2 \ + \n[\ + I^D). 

Let us first deduce that the operator R = AM A' 1 is defined on the domain C A - If 

10) e C Aj 



Dom (i?) ^ Mi^ x |0) G D A ^ ^ | (in'lMi^ 1 ! 



< oo 



E 



^{An'\M\A~ l n)(n\ 



< oo «=> 



E 



< oo. 
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Let us show now the last estimate is actually fulfilled. If \<f>) G L 2 (T 2 ^ 



y]Rn> tn (n\ 



< 



^ T,(j2\ R <n\\(n\4>)\) 

n' \ n / 

||0f ^ e - 2 <KI+KI) (j2 e ~ c{lni 

n' \ n 



l+lnal) 



< OO. 



Therefore the operator -R = AM A 1 is defined on the domain Ca and its matrix elements 
are obviously (n'\R\n) = R n i >n . With the same arguments as those used earlier for the 
expanding map on S 1 , we deduce that R extends to a trace class operator on L 2 (T 2 ). 

4.2.1 Proof of the exponential concentration of R.P. resonances 

The proof of Theorem 7 is very similar to the proof we gave in Section 3.3. Here we use 
the same notations and emphasize the differences. 



Let R be the trace class operator obtained in Theorem 6, with the estimation 
exp (— c(|ni| + | rz-2 1 + \ n [\ + l n 2l)) ° n its matrix elements, with c > 0. 



(n'\R\n) 



< 



Lemma 2. Let fij, j = 0, 1, . . . be the non zero singular values of R, such that fij + i < fij, 
repeated as many times as the value of their multiplicity. Then 



l'.i 



(23) 



with C\ > 0. 

Proof. From the min-max theorem, 



u ? = mm max 

VcL 2 (T 2 ),codimV=j veV, \\v\\=l 



Rv 



Consider the Fourier basis \n),n = (ni,n 2 ) G Z 2 and Vi = span(|n)) 
codimV/ = (2/ + l) 2 . If \v) G V/, we compute 



max(|ni|,|n 2 |)>Z' 



hence 



\n')(n'\R\n)(n\v) 

n'&L\n\>l 

< \H ( ex p(- c (l n 'il + K2D) J 

Vn'eZ 2 / 

< \\v\\Ce- d , C>0 

We deduce that ^ < C x e~ c ^l 2 , with C x > 0. 



Y exp(-c(|m| + |n 2 |)) 

, n/ max(|ni|,|n2|)>i 



□ 
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For non zero eigenvalues, define 

lj = \og\\j\, mj = log iij 

(These sequences tend to — oo as j —> oo). Inequality (19) reads Y^j=oh — X]J=o m i- ^q. 
(23) gives rrij < logCi - cVJ/2. We deduce that j™=oh < (n + l)logCi - f E"=o v^- 
But E;=o v^J > Jo" v 7 ^ = §n 3 / 2 . Hence ^ E;=o ^' < kgd " f (SSy- But Z n < /, for 
j < n, so /„ < Ej=o^ ^ logCi - f (i+i/n) ' which P roves Theorem 7. 

5 Conclusion 

For specific models of chaotic dynamics, namely real analytic expanding maps on the circle 
S 1 and real analytic hyperbolic maps on the torus T 2 , we have shown that the decay of time 
correlation functions can be described by a trace class operator in L 2 (S 1 ) (resp. L 2 (T 2 )). 
We have followed an approach similar to the "complex scaling method", well-known to 
study the decay of quantum states in open quantum systems. As explained in the intro- 
duction, this approach has been already used by V. Baladi and M Tsujii [4] for hyperbolic 
diffeomorphisms in a more general context, but our methods differ sightly and allowed us 
to obtain different results. To make a more precise comparison, our operator A defined in 
Eq.(ll) and the conjugation in Eq.(12), correspond to their definition of anisotropic norms. 
But they use powers of the Fourier components whereas we use exponential of them. This 
exponential is important for us to obtain a trace class operator R, as explained in section 
3.2.1. 

In this paper, some "semi-classical aspects" of hyperbolic dynamics have appeared many 
times: (i) in Theorem 10 concerning the localization of the matrix elements in Fourier 
space, and the remark which follows, (ii) in a remark on the semi-classical Weyl law after 
equation (13), and (iii) with Figure 5 where a cluster of eigenvalues suggests some semi- 
classical tunnelling effect. The role of semi-classical parameter is played by the inverse of 
the distance in Fourier space: H = 1/ \n\. A direction of research would be to make this 
semi-classical theory more precise. 

We would like to comment some limitations of our results and possible extensions of 
them. First we have assumed that the dynamics is given by a real analytic map. This is 
a severe limitation because in hyperbolic dynamical system theory one has to use Holder 
potential functions [3] [14]. We have presented here the results in their simplest form. In 
particular, for expanding maps on the torus, we have shown that R is trace class for any 
expanding map, but for hyperbolic maps on the torus we have assumed that the non linear 
perturbation / is weak enough. In geometric terms we have supposed that is weak 

enough so that the unstable and stable foliations are respectively contained in fixed cones 
(the cones adapted to the linear map, and which enter in the definition of A eq.(ll)). 
It could be possible to generalize in this direction and treat in this way any uniform 
analytic hyperbolic map on the torus, using a local choice of cones. As in [4], this could be 
possible using pseudo-differential operators instead of A. Then the localization property 
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of the matrix elements in Theorem 10 would be replaced by a "microlocal" version in the 
cotangent space T*S 1 (iesp. T*T 2 ). Some other directions of research could be to take 
advantage of the relative simplicity of this approach to investigate non uniform hyperbolic 
dynamics or other kinds of dynamical systems which exhibit some chaotic behaviour, where 
many questions remain open. 
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